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Well known from the sixties, the pressure of e. g. massless 
theory may be written as a series of 2PI-diagrams (skele¬ 
tons) with the lines fully dressed. Varying the self-energy 
n in this expression, it turns into a functional U \Y] having 
a maximum in function space at V = 77. There is also the 
Feynman-Jensen thermal variational principle a poten¬ 

tially non-perturbative tool. Here actions S are varied. 

It is shown, through a few formal but exact steps, that the 
functional U is covered by V. The corresponding special 
subset of trial actions is made explicit. 


I. INTRODUCTION 

Beyond perturbation theory, we are in search for some 
optimization calculus. The thermal variational principle, 
while well appreciated in non-relativistic quantum statis¬ 
tics, plays a merely minor role in field theory so far. At 
present, there are two apparently different such princi¬ 
ples, which we shall call the Luttinger-Ward principle |^] 
and the Feynman-Jensen principle |^,||. For the study 
of the possible difference of these two principles we shall 
concentrate on scalar (f>‘^ theory in the first two sections. 

In 1960 Luttinger and Ward Q made a remark in 
parentheses, that the pressure p of their fermionic sys¬ 
tem, after the diagram lines were fully dressed, be¬ 
comes minimal under variation of the self-energy. For 
the massless thermal (()^-system with Lagrangian C = 
(dcj))'^/2 — this property may be recapitulated 

as follows. The skeleton version for the pressure p, which 
is —1/V times the free energy F, reads 

p = Ayi^(nip)-p^) + ly ^ 

y 2y \ ^ > ’ 2y n{p) - p^ 

I oo 

+ r with rsjfpr ■ (1) 

^ n—1 

n{P) is the exact self-energy, J2p = I'^^P 

is the thermal sum-integral, and fy are the 2-particle 
irreducible contributions of order 5^" to the logarithm 
In(Zint) of the partition function: 


r = 3 OO +12 dD + 2 • 12^ ( 2 ) 

+ ^-123 QQ -b 6-123 QQ +... ^ (2) 

The result (0) is found in J0j^]. It can be derived by a 
Legendre transformation |]6[0 or, equivalently, by mini¬ 
mizing the free energy ||]. It was recently taken up in , 
see also the article by A. Peshier in these proceedings. 

Three minor modifications make the above expression 
(|^) to become the functional U of what we call the 
Luttinger-Ward variational principle. We multiply (|^) 
with PV . We supply the volume V with periodic bound¬ 
ary conditions, hence the thermal sum turns into the 
’’bare sum” Op — On Op = Op • Third, we replace 
the self-energy 7T (P) by some function Y (P) to be var¬ 
ied and introduce the notation G{P) = 1/{Y{P) — . 

Then, as is seen shortly, the principle states that 

U[Y] = ^^ In (G(P)) + i £ P(P)G(P) 

P P 

+ r[G] < pvp . ( 3 ) 

The fact that here the pressure is at maximum, in con¬ 
trast to that of [Q, is due to the boson content of our 
system. The replacement P —> V had to be performed 
even in the last term, i.e. in the diagram lines of (0). 
So, these lines, the bones of the skeleton, have become 
variable propagators G(P). 

For taking the functional derivative of U with respect 
to Y{Q) at an arbitrary ’’position” Q, we appreciate the 
general functional relation 0 

2GlSGof = G (/ = ln(Z)) , (4) 

between bare and dressed lines, which is valid separately 
in each loop order |0|. Hence 

2PI 

2G‘^6gT[G] =G2P^[F] = 

= _G2p2PI[y] (5) 

and 


2 Sy^Q) U = -G^ (y - p2PI[y]) 1 0 
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because the full n{Q) (up to a given order n) is recon¬ 
structed by iterative use of Ym = with m < n 

1^. At any position Q, when Y increases there, Sy{q)U 
reaches the zero coming from positive values and turning 
into negative. So, the functional U has a true maximum 
at r = 7T. 

Some warning is in order, if the sum ^ in U is 
truncated. Let 2n' be the highest g-power contained in 
U . Then, at best, the resulting 77 would be correct to 
order only (n' >4). 77 in for instance, needs 

U to be developed to This unfortunate fact 

is due to reduction of 5 -order by evaluation. Things are 
realized by power counting of soft scale contributions, or 
along the lines given by Braaten and Nieto [^. One 
might question for a better functional, which is free from 
this defect. For the following we stay with the exact 
expression (^. 

To introduce the Feynman-Jensen variational principle 
we must distinguish between the action S' of a trial theory 
and the action S, of the true (()^-system under study. 
With or without a bullet index, we have by definition 
S = -fc,f = J^dT Jd^r and f3=l/T. Then g 




J V(j) e-^' = (^j V(j) 


JV(j)e 
f Vfj) e~^ 

(7) 


with Z, Z, the partition functions of trial and studied 
theory, respectively. The inequality in is the famil¬ 
iar Jensen inequality jig applied to functional integrals. 
The functional measure is obviously irrelevant here. Note 
that the average refers to the trial action S, which usu¬ 
ally might be the simpler theory. Taking the logarithm 
of ( 0 ) we arrive at 

7;[S] = ln(Z) + (S-S.) < ln(Z.) . ( 8 ) 


system with some momentum-dependent mass Gq^ + 
such that a given exact propagator G = 1/{Y — P^) turns 
out perturbatively. Note that 77 = G~^ — Gq^ differs 
from Y = G~^ + P^ . Since F[G] is a functional of G 
only, all the bare lines, which constitute the partition 
function Z diagrammatically, must be re-expressed by 
G: 




n—l 


\n(Z)-^^\n{G)-^^nG ^ 


G 


ln(Z*)-^£ln(G)-^E(l 


■ (9) 


Without star index, the inner line of (g) is nothing but 
the skeleton formula (g, but this time for the mas¬ 
sive system just introduced. The prescription has the 
following definite meaning: solve the relation d), i.e. 
2 Gq(5go/ = G , for Go ; then replace all Go’s in the square 
bracket by the resulting functional Go [G] = Go * . This 
is a Legendre transformation |^g|. As G does not change 
under the ^-operation, the second term reached the sec¬ 
ond line unchanged. But it will now be canceled when 
using (g in (g : 

U[Y] = ^^YG +In 

< PVp = ln(Z.) . (10) 


The inequality refers to variation of Y. While still G = 
l/(y — P^), the objects Go* and Z* are some nontrivial 
functionals of G, hence of Y . 

Step two. We transcribe the inequality into func¬ 
tional integral language. Exponentiating we have 


Also the functional V has a true maximum, this time 
at S = S, . The principle (g was studied with detail in 
g and applied to gauge fields. Working with free trial 
theories, several known results could be reproduced in 
g. But to get more information from the principle, the 
free trial space has turned out being too poor. 


II. THE WAY U IS CONTAINED IN V 

In three steps the principle U will be reformulated 
until it has the structure of V . 

Step one. To get rid of the inhnite sum T in i), 
we exploit the fact that the l.h.s. of g) regains a 
physical meaning by the replacements YG —> PG and 
U ^ \n (Z). Here, Z is the partition function of a 


Lf'(YG-l+^) Y ^ Y 

gi£(^G-i-i-^) < Jv(j,e~^- , (11) 

where S, = -/^(£o + Ant), but S'* = -f^(£o* + 

+£int) with 

-^0 -^o|Go—*G o* ; 

- = = So* . ( 12 ) 

We turn to the exponential in ( 0 ), which might be 
some functional average over two fields cj). In fact, if we 
first separate a factor G by 
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LjO* Vt-fO* ^00 


where Goo = and write this factor as 

G(F)2G^o*SGoJn(Z^) = 


^ <5go* I 


then the exponential becomes 


(13) 


(14) 




1 

Go* 


1 

Goo 


)^»). 


= {So,:-So,)^ = iS^-S.' 


(15) 


Using ( p^ in (0 , the result 

< Z, = (16) 

formally agrees with (^ and is a Jensen inequality. 

Step three is merely one in mind rather than in for¬ 
mulation. So far we varied functions Y around 11 or, 
equivalently, G’s around 1/(77 —P^). But this is not the 
variation in the Feynman-Jensen functional. The latter 
depends on actions. So, we now change the philosophy 
and consider Gq* to be the variable function. Hence 
& becomes the Feynman-Jensen principle, indeed. But 
the class of trial actions is very restricted. The full 
interaction Ci^t is part of and remains untouched. 
Due to (0 S'hc differs from S', only in some variable 
momentum-dependent mass term. Moreover, the solu¬ 
tion to the optimization problem has become trivial: the 
equality sign in (16) is simply reached at vanishing mass 
term. The true value of using V, in contrast, is in opti¬ 
mizing approximations for a physics too hard for an exact 
or even perturbative analysis. 

From the fact, that we are on a very safe ground with 
the special Feynman-Jensen variational principle dl), 
we are led to reverse the order of the above three steps. 
Then they guide a possible derivation of the skeleton for¬ 
mula (|^ including its variational property (^) . This idea 
is followed up in the next section. 


fixing term, a term arising from the Faddeev-Popov de¬ 
terminant and the 3- and 4-point interactions, 


S, = So + S^.i. + SFp[A] + Si^t . 


(17) 


For the functional inequality ( [Tq ) to make sense (but of 
no relevance in the sequel), 5 'fp must be viewed as a 
functional of the fields [|j. The trial action Sit: and 
the action S', have the same gauge fixing parameter a 
(see § II.B of Q). Then, the only difference between the 
two actions is a momentum-dependent mass term, which 
we immediately write as 


S*-S, = 




Go = 


1 


AliP) 


(18) 


The fully dressed Greens function of the ^-theory 
enters when averaging (0) 




= ?/^^E 


LGo*(P) Go(P)J 


Gt.u[P) ( 19 ) 


with n = Ga = — 1 ■ Introducing a variable matrix 

self-energy Y^'^(P) we are led to 


1 


G, 


0 * 


1 


G = 


1 1 
Go * Go 
1 


1 


Go"' - 


G, 


G-l-YG 


0 * 


( 20 ) 


with all products being Lorentz matrix multiplications. 
Some signs differ compared to (0 due to Minkowski 

metrics H-and the conventions of |3j. Chan ging 

from variable Go it: to variable Y, the logarithm of ( |lq ) 
reads 


U[Y]=n^X[G] - 5£uG + r[G 


< ln(^,) 


( 21 ) 


with 


r[G] = 


\n{Z)-n^X[G] + ^^ 


2 V VGo 


G- 1 


III. THE SKELETONS OE YANG AND MILLS 


= E/> 


2PI 


( 22 ) 


We derive the skeleton formula for the pure gluon sys¬ 
tem by following the above equations in backward di¬ 
rection. In covariant gauges the action includes a gauge 
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So far, the term with Ar[G] has been only added and 
subtracted. But now we require the square bracket to 
be the sum over 2PI contributions to In (Zsi:). With the 





















Yang-Mills counterparts of (i) and (|), as listed in ([2^, 
this condition turns into the matrix differential equation 
in the third line : 


2Go{SGof))Go = nG , 

Sg^X[G] = ^ . (23) 


The solution to the X equation needs (^) and is 
given (apart from some Y-independent constant) in ( p^ ) 
below: first term. Herewith, one might drop the question 
mark in (|^). 

To supply the first two terms of ( 0 ) with detail, we 
need the structure of the exact gluon propagator. Due 
to Weldon [Q this propagator and the self energy may 
be written as 

yz 

_ B'"") , (24) 


where At^e = 1/(T’^ — Yt^i). Apart from A = g — B — D 
all the Lorentz matrices in (p4|) are dyadic products of 
the four-vectors P = {Po,p) and P = {p, PqP/ p), or 
oi R = P — yP and R = P — yP , respectively. To 
be specific: B = —P o P/P'^ , D = P o P/P"^ , B = 
—R o R/P"^ , D = Ro R/P^ . We emphasize that even 
the Weldon coefficient y is varied. With (^^, and after 
preparing the trace over all pairs of the above matrices, 
one verifies = 2YtAt + Y^A^ — ay^ with ease. 

The Luttinger-Ward functional U for the gluon system 
then finally reads 


U[Y] = ^^\n 

p 


1 


[Yt-P^l 

Ye-P^ 


n / 21) 

2 ^ [Yt-P^ 


Ye 

Ye-P^ 



+ r 


(25) 


The maximum is reached at Y( = iTj, 1) = 77^ and y = b, 
where the coefficient b is defined by (|^) at Y^"^ = 
and is given by \/2 lid/He in earlier notations jlj]. At 
maximum ^ U becomes the pressure of the hot gluon 
medium. 


of trial actions is reduced to a rather trivial variation of 
mass terms. Nevertheless, the machinery relating the two 
principles can be used to derive the skeleton formula for 
the pressure of both, the massless scalar (f>‘^ theory and 
the gluon system. 

I thank Jens Reinbach for deriving the Luttinger-Ward 
functional for Yang-Mills fields along the lines given 
in§. 
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IV. CONCLUSIONS 


The competition between the Luttinger-Ward and the 
Feynman-Jensen variational principle is won by the lat¬ 
ter. It contains the former as a special case, whose space 
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